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Abstract 

> . 
Tilj- , The improved calculation method of quark distributions in hadrons in the frame- 

^N I work of QCD sum rules is presented. The imaginary part of the virtual photon scat- 

f— ^ ' tering amplitude on some hadronic current is considered in the case, when initial and 

final virtualities of the current pf and p^ ^^^ different, pf 7^ pi- The double Borel 
^^ I transformation in Pi,P2 is applied to the sum rule, killing background non-diagonal 

transition terms, which deteriorated the accuracy in the previous calculations. Va- 
lence quark distributions in the pion were found in a good agreement with the data, 
determined from the Drell-Yan process. 



In my talk I present the recent results in determination of quark distributions in mesons 
by using the improved QCD sum rule approach, which were obtained by A.Oganesian and 
myself |]I|. I recall the idea of the method 0, ^. Consider the 4-point correlator of two 
electromagnetic currents and two currents with quantum numbers of some hadron (for 
clarity the axial current, corresponding to charged pions is considered): 

(1) 

J5A = u-f5-fxd (2) 

where qi,q2 and Pi,P2 are the photon and hadron current momenta, correspondingly. 
Assume, that g^ and p"^ are negative, qf = q2 = 1^, \ l"^ \^\ Pi \y \ P2 \ ^^^ \ Pi l~l P2 1^ 
R~^, where Re is the confinement radius. As was shown (|^), in this case the imaginary 
part in s-channel of the correlator (|I|) is dominated by small distances in all channels. 
Particularly, at Pi = P2, qi = Q'2 the closest to zero singularity of Im Hf^^Xp in the t-channel 
is given by {x = — g^/2z/, u = pq) 

t = -4-^p^ (3) 
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and t ~| p^ 1^ Rc^ at nonsmall x, unlike the case of Ilfj,u\p itself (the forward scattering 
amplitude), where t = (for massless quarks) and large distances in the t-channel are 
of importance. This fact allows one to use operator product expansion (OPE) in 1/p^ 
for calculation of Im n^j^^p in QCD, besides the standard OPE in 1/q^. As follows from 
(^ the approach does not work at small x. It does not work also at large x, close to 1. 
Physically it is evident, because this is the resonance region. 

Phenomenologically Im Il^yXp is represented by contributions of physical states using 
the double dispersion relation in p^, p|. Among various tensor structures of li^uXp it is 
convenient to consider the structure {P^PyP\Pp/u)Ii{j)\,p\,q^,x)^ where P = {pi+p2)/2, 
X = —q^/2{qiP). The double dispersion representation of Im 11 has the form: 



r fr/ 2 2 ^ / ^ , ffi^^u) [ V{x,u) f f p{x,Ui,U2) ,. 

Im ll(pi , po, xj = aixj + / 5-dM + / 5- + / dui / 0^27 ^77 ^ (4) 

{ u-pi I u~pi J 7 {ui - pi){u2 - pi) 

(the q^ dependence is omitted). In the previous treatment of the problem the dis- 
persion representation (H) was considered in the limit pi = pi, = p^ and the single Borel 
transformation in p^ was performed. This results in appearance of nondesirable back- 
ground terms, which deteriorate the accuracy of calculations and in the case pion even 
do not allow one to find the quark distributions at all. The spectral functions in (^ were 
represented by contributions of the lowest state (pion) and continuum 



p(mi, M2, x) = fl- 2-kF2{x)6{ui - ml)6{u2 - ml) + p^{x)0{ui - so)9{u2 - sq) 

ip{x, u) = ipi{x)6{u - ml) + ip2{x)9{u - sq) (5) 

where /^r = 131 MeV and F2{x) is the pion structure function, sq is continuum threshold. 
(pi{x) corresponds to the contribution of the nondiagonal transition 7* + vr ^ tt* + 7* and 
is unknown. It is not suppressed in comparison with the main term ~ F2{x) by the single 
Borel transformation and must be accounted in the final sum rule. 

The lowest order term of OPE corresponds to the box diagram of Fig.l. At pi = p2 
its contribution to Im, U is equal 

3 1 

/mn(p^,x) = rx^(l — x) (6) 

7rp^ 

Using (^), (|^) and (^) and performing the single Borel transformation in p"^ we get the 
sum rule (M^ is the Borel parameter) 

^x2(l - x)(l - e-«/^^^) = 27r/>.(x)^ + ^i(x), (7) 

where Ut^{x) is the distribution of valence u quarks in the pion (the pion mass is neglected). 
Looking at M^ dependence in (0) it becomes evident, that in this approach the attempt 
to separate the pion contribution from the background by studying M^ dependence (e.g. 
differentiation over 1/M^) is useless - up to small correction ~ g-so/M ^j^g ^^^^ diagram 
contributes to the background only. 



Let us now improve the method by considering the case p^ 7^ p\ and performing 
the double Borel transformation in Pi,P2- The double Borel transformation kills the 
nondesirable first three terms in (^. Instead of @ we have now 
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n(p?,P2)3;) = —x^(l — x) du du'- 7T-—^ ^ (8) 

TT ' 'II {u-pI){u'-pI) 

and after double Borel transformation with parameters M^, M2 the sum rule arises 



uA^) = ^,^x{l - x){l - e-^^^/'% (9) 



3 M2 
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where it was put Mf = M2 = 2M'^. The calculation of the pion decay constant /^r, 
performed in H in the same approximation leads to 



/.' = ^M2(l-e--/M^) (10) 

The substitution of ( [T0| ) into (H) gives 

■u^(x) = 6x(l — x) (11) 

Therefore the necessary conditions 

1 

UT^{x)dx = 1, / xuT^{x)dx = 1/2 (12) 


are fulfilled in this approximation. 

Account now the perturbative and nonperturbative corrections. Restrict ourselves to 
LO perturbative corrections, proportional to ln{Q'^/M'^) and choose Q^ = Qq = 2 GeV^ 
as a point at which we calculate the quark distributions. The LO perturbative correction 
results in multiplication of bare loop distribution (|]) by the factor 



r{x,Q' 



as(Myn(QyM^) ^ , „ , , 2(l-2a;)/ra. 
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(13) 



The higher order terms of OPE in l/p^ starts from contribution of gluonic condensate 
(0 I ^ui/^ui/ I 0) of dimension 4. The calcukation was performed in the Fock-Schwinger 
gauge x^A^{x) = using the program of analytical calculation REDUCE. Surprisingly, 
the sum of all diagrams, proportional to gluonic condensate vanishes after double Bore- 
lization and gluonic condensate does not contribute to the sum rule. 

One-particle irreducible diagrams, (no loop diagrams) resulting in appearance of 
(5(1 — x) in Im U and the diagrams, arising from their QCD evolution are disregarded 
because the calculation method is inapplicable at x = 1. There are two vacuum expecta- 
tion values of dimension 6: {g^ f"'^'^G'^^j^G\,yG')^^ and (V'V')^ (for V'O-?/' ■ V'O?/' operators the 
factorization hypothesis was used). The calculation shows, that the contribution of G^ 
operators vanishes after double Borel transformation and only {ijjipy contribute at d = 6. 
The fonal result for valence w-quark distribution in vr^ with account of LO perturbative 
corrections and OPE up to (i = 6 is given by 



3 M2 2,, , 
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'^^^^g2^.(^_g-so/M^)_47ra,(g2) .47ra,(M2)a2 uj{x) 



(27r)4 ■ 37 • 26 . M6 x3(l - x)3 

(14 
where uj{x) is the polynomial of 4-order in x, a = — (27t)'^ (uu) . The analysis of the sum rule 
(14) shows, that it is fulfilled at 0.15 < x < 0.7; the power corrections are less than 30% 



and the continuum contribution is less than 25%. Stability in M^ at 0.4 < M^ < 0.6 GeV^ 
is good. The final result for xu.„{x) (at M^ = 0.45 GeV^ and Sq = 0.8 GeV"^) is shown in 
Fig. 2. Fig.2 shows also the curve of the valence u-quark distribution in the pion, found 
in from the fit to the Drell-Yan procces data. The agreement is satisfactory, especially 
bearing in mind, that nonaccounted NLO perturbative correction would increase Uj^lx) at 
small X and decrease at large x. If we make the assumption that at x < 0.15 UTr{x) ~ 1/y/x 
according to Regge behavior, and at x > 0.7 M7r(x) ~ (1 — x)^, then the moments can be 
found 

1 1 

Mo = / u^{x)dx ^ 0.84 Ml = / xun{x)dx ^ 0.21 (15) 
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Figure 1: Diagrams, corresponding to the unit operator contribution. Dashed hues with 
arrows correspond to the photon, thick sohd - to hadron current 
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Figure 2: Quark distribution function in pion, noted "total". For comparison fit from 
[5], noted "GRV", is shown. Also bare loop ("bare") and bare loop with nonperturbative 
corrections (noted "1"), are shown 



